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Abstract 

We study the ground-state (gs) properties of the frustrated spm-i Ji~J2-J?, Heisenberg model on 
the two-dimensional honeycomb lattice with ferromagnetic nearest-neighbor (Ji = — 1) exchange 
and frustrating antiferromagnetic next-nearest-neighbor (J2 > 0) and next-next-nearest-neighbor 
(J3 > 0) exchanges, for the case J3 = J2. We use the coupled-cluster method implemented to 
high orders of approximation, complemented by the Lanczos exact diagonalization of a large finite 
lattice with 32 sites, in order to calculate the gs energy, magnetic order parameter, and spin-spin 
correlation functions. In one scenario we find a quantum phase transition point between regions 
characterized by ferromagnetic order and a form of antiferromagnetic ("striped") collinear order 
at J2 ~ 0.1095 lb 0.0005, which is below the corresponding hypothetical transition point at J2 = y 
(~ 0.143) for the classical version of the model, in which we momentarily ignore the intervening 
noncollinear spiral phase in the region < J2 < \- Hence we see that quantum fiuctuations appear 
to stabilize somewhat the collinear antiferromagnetic order in preference to the ferromagnetic order 
in this model. We compare results for the present ferromagnetic case (with Ji = —1) to previous 
results for the corresponding antiferromagnetic case (with Ji = +1). The magnetic order parameter 
is found to behave similarly for the ferromagnetic and the antiferromagnetic models for large 
values of the frustration parameter J2. However, there are considerable differences in the behavior 
of the order parameters for the two models for J2/\Ji\ ^ 0.6. For example, the quasiclassical 
collinear magnetic long-range order for the antiferromagnetic model (with Ji = -|-1) breaks down 
at ~ 0.60, whereas the "equivalent" point for the ferromagnetic model (with Ji = —1) occurs 
at J| ~ 0.11. Unlike in the antiferromagnetic model (with Ji = +1), where a plaquette valence- 
bond crystal phase intrudes between the two corresponding quasiclassical antiferrmagnetic phases 
(with Neel and striped order) for < J2 < J2^, with ~ 0.47, we find no clear indications 
at all in the ferromagnetic model for an intermediate magnetically disordered phase between the 
corresponding phases exhibiting ferromagnetic and striped order. Instead the evidence for the 
ferromagnetic model (with Ji = —1) points to one of two scenarios: either there is a direct first- 
order transition between the two magnetically ordered phases, as mentioned above; or there exists 
an intervening phase between them in the very narrow range 0.10 ^ J2 ^ 0.12, which is probably 
a remnant of the spiral phase that exists in the classical counterpart of the model over the larger 
range < J2 <\. 
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I. INTRODUCTION 



In recent years frustrated quantum spin systems on regular two-dimensional (2D) lattices 
have aroused a great deal of research interest.-"- In particular the interplay of magnetic frus- 
tration and quantum fluctuations has been seen to be a very effective route to destabilize 
or destroy magnetic order and thereby to create new quantum phases. Such 2D magnetic 
systems can thus in turn develop a diverse array of phases with widely different ordering 
properties, such as antiferromagnets with quasiclassical Neel ordering, quantum "spirals", 
valence-bond crystals/solids, phases with nematic ordering, and spin liquids. Other factors 
that influence the ground-state (gs) phase structures are the nature of the underlying crystal- 
lographic lattice, the number and nature of the bonds on this lattice, and the spin quantum 
numbers of the atoms localized to the sites on the lattice. The theoretical investigation of 
these models has proceeded hand in hand with the discovery and experimental investigation 
of ever more quasi-2D magnetic materials with novel properties. 

One of the most intensively studied of all of the frustrated 2D models is the spin-i 
J1-J2 Heisenberg antiferromagnet (HAF) on the square lattice with nearest-neighbor (NN) 
bonds (of strength Ji > 0) competing with next-nearest-neighbor (NNN) bonds (of strength 
J2 = aJi > 0). This quantum system has two different quasiclassical phases with collinear 
magnetic long-range order (LRO) at small {a < ~ 0.4) and large {a > a^^ ~ 0.6) values 
of the frustration strength parameter a, separated by an intermediate quantum paramagnetic 
phase with no magnetic LRO in the regime < a < ac2- Interest in this model has 
been greatly stimulated recently by its experimental realization in such layered magnetic 
materials as Li2VOSi04,-'^ Li2VOGe04,- and VOM0O4.- The syntheses of such layered 
s has stirred up a great deal of renewed interest in the model (and see 



quasi-2D materia 

also, e.g., Refs. 7H10|). Amongst several methods that have been very successfully applied 
to the J1-J2 model has been the coupled cluster method (CCM),— "— which has also been 
applied to many similar strongly-interacting and highly frustrated spin-lattice models with 
comparable success. Other frustrated 2D models that have similarly engendered great recent 
interest include the spin-i HAFs on the triangular—*^ and kagome lattices.—*^ 

There has been a large amount of recent experimental investigation of the properties 
of quasi-2D magnetic materials with a ferromagnetic (FM) NN coupling (Ji < 0) and an 
antiferromagnetic (AFM) NNN coupling (J2 > 0). Examples include Pb2VO(P04)2,— "— 
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(CuCl)LaNb207,-^ SrZnVO(P04)2,-'2^'26 BaCdVO(P04)2,— '^^^ PbZnVO(P04)2,^ and 
(CuBr)LaNb207.— These experimental studies have also served to reignite interest in the 
theoretical investigation of the gs and thermodynamic properties of the FM J1-J2 model, i.e., 
the model with FM NN exchange ( Ji < 0) and frustrating AFM NNN exchange ( J2 > 0).-^ — 
Interestingly, arguments for the existence of a spin-nematic phase between two quasiclassi- 
cal magnetically-ordered phases were presented.— i^i^i^ On the other hand, the existence 



of such a non-classical magnetically-disordered phase was also questioned in Ref. 



39|. 



Other systems that have grown in importance in the last few years are various spin- 
I magnetic models defined on the 2D honeycomb lattice. Several such systems have been 
both theoretically and experimentally studied^""— intensively, partly because of their special 
properties and partly due to the recent syntheses of various quasi-2D honeycomb-lattice 
materials. One reason for the theoretical interest in such models on the 2D honeycomb 
lattice is that a spin-liquid phase has been found for the exactly solvable Kitaev model,— in 
which the spin-| particles reside on just such a lattice. Furthermore, the honeycomb lattice 
is obviously germane to the very active research field of graphene, where the relevant physics 
may well be described by Hubbard-like models on this lattice.— "^"^ Interestingly, Meng et 
alr^ found that for the Hubbard model on the honeycomb lattice with moderate values 
of the Coulomb repulsion parameter f/, strong quantum fluctuations lead to an insulating 
spin-liquid phase between the non-magnetic metallic phase and the AFM Mott insulator 
phase. From the experimental side recent observations on the spin-| honeycomb-lattice 
HAF Bi3Mn40i2(N03) demonstrate a spin- liquid-like behavior at temperatures much lower 
than the Curie- Weiss temperature.— 

We have recently studied^ the AFM honeycomb model for the case where the 

spin quantum number s of each of the spins on every lattice site is s = |, and with AFM 
nearest-neighbor exchange bonds (Ji > 0) in the presence of frustration caused by AFM 
NNN bonds (J2 > 0) and with next-next-nearest neighbor (NNNN) bonds of strength J3 
also present, for the special case where J3 = J2- We found^ that the scenario of deconfined 
criticality may hold for this model (and see also Ref. 48|). To date there exist only limited 
studies of the corresponding FM model (namely where Ji < 0). In this paper, 

we further the investigation into the FM honeycomb model with FM NN bonds 

(of strength Ji < 0) in the presence of frustrating AFM NNN bonds (of strength J2 > 
0) and NNNN bonds (of strength J3 > 0). Once again we consider only the interesting 
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special case where J3 = J2. We focus our attention in the present study particularly on 
the detection and characterization of the gs phases of the quantum model. Bearing in 
mind the controversial discussion of the corresponding J1-J2 square-lattice model with FM 
NN exchange bonds (Ji < 0), the question naturally arises as to whether any indications 
for a non-classical magnetically-disordered phase might now be found for the honeycomb 
model. To determine the relevant gs phases and their properties we calculate the gs energy, 
the spin-spin correlation function, and the magnetic order parameter for the stripe-ordered 
state discussed below that is present as a gs phase in the corresponding classical version 
(equivalent to the s — j- 00 limit) of the model. 

In view of its proven past ability to give results of high accuracy for a wide variety of 
highly frustrated 2D spin-lattice models, we again use the coupled cluster method (CCM) as 
our main computational tool in this paper. Additionally, we use the exact diagonalization 
(ED) method for a large finite lattice of = 32 spins as a validity check of our CCM results. 
Since at the classical level the model now under consideration also exhibits some similarities 
with the corresponding model with AFM NN bonds (Ji > 0), we compare our results for 
the quantum model of the FM case ( Ji = —1) with those of the corresponding AFM case 



The rest of the paper is organized as follows. After describing the model in Sec. [TTl we 
apply the CCM to investigate its gs properties. The CCM itself is very briefly described in 
Sec. mil before presenting and discussing our CCM and ED results in Sec. HVl We conclude 
in Sec. |V]with a summary of the main results. 

II. THE MODEL 

The Hamiltonian of the spin-| Ji-J2^Jz Heisenberg model on the honeycomb lattice, 
which we studied recently^ for the AFM case (Ji > 0) is defined as 



where i runs over all lattice sites on the lattice, and where j runs over all NN sites connected 
to site i by Ji bonds, k runs over all NNN sites connected to site i by J2 bonds, and / runs 
over all NNNN sites connected to site i by J3 bonds, but counting each bond once and once 
only in the three sums. Each site i of the lattice carries a spin-| particle with spin operator 




(1) 



6 













1 









(a) FM 



(b) striped AFM 



FIG. 1: (Color online) The J1-J2-J3 honeycomb model, showing (a) the ferromagnetic 
(FM) state and (b) the classical striped antiferromagnetic (AFM) state. The arrows 

represent spins located on lattice sites •. 



Sj = (sf,sf,sf). We note that precisely the same model has also been studied recently on 
the square lattice, both in the case where all the bonds are AFM in nature,— and in the FM 
case where Ji < and J2 > 0, J3 > 0.— 

The aim of the present work is now to study further the spin-^ Ji-J2^J3 FM model 
(namely the above model in the case Ji < 0) on the honeycomb lattice.— "—1^ The lattice 
and the exchange bonds are illustrated in Fig. [U^a). The classical gs phase diagram for the 
J1-J2-J3 AFM model (with Ji > 0) on the honeycomb lattice model displays coUinear Neel 
and striped phases, both AFM in nature, as well as a spiral phase. These phases meet in a 
triple point at J3 = J2 = Ji/2 (and for more details see, e.g., Ref. 45|). For the remainder of 
this paper we again focus on the case where J3 = J2 > 0, but now where the NN exchange 
bond is FM in nature (Ji < 0). The gs energies of the only two corresponding classical 
collinear states are then given by 



N 

striped AFM 

N 



3 X 9 ^ 
- Ji + -J2 
2 2 



(2) 



for the FM state and collinear striped AFM state shown in Figs, [^a) and (b) respectively. If 
these were the only gs phases in this Ji < regime we would thus have a classical transition 



between the FM state and the striped AFM state at = -yJi 0.143 for Ji = -1) 
and a classical energy per site at this point of E^^/N = — ^ ^ —0.214 for the s = | system 
with Ji = —1. For the corresponding AFM model with Ji > such a striped AFM state 
also exists as stated above, but the classical transition between the AFM Neel state and the 
striped AFM state is at J2 = The reason why the corresponding phase transition in 



the FM model occurs at a smaller value of the frustration parameter, J^^Vl^^il — y than the 
value J2/J1 = I for the AFM model is due to the J3(> 0) NNNN exchange bonds that act 
to frustrate the fully polarized FM state, whereas they reinforce the AFM Neel state. By 
contrast, the J2(> 0) NNN exchange bonds act to frustrate the Ji bonds for both the FM 
state of the FM model and the Neel state of the AFM model. We note that the classical FM 
state is also an eigenstate of the Hamiltonian, with energy eigenvalue equal to the energy of 
its classical FM counterpart. 

We note, however, that in fact the classical Heisenberg model on the honeycomb 

lattice with Ji < 0, J2 > 0, and J3 > also has a spiral phase that intervenes in a very 



narrow strip between the FM phase and the collinear striped AFM phase. (In Refs. j45 
this is referred to as phase V.) The region in the x-y plane (where x = J2/J1 and y = J3/J1] 
in which it is the stable gs phase in the case Ji < is bounded by the three curves (i) y = 



— \ < X < — |, (ii) y = — |x — |, — I < X < |, and (iii) y = ^[1 — 6x — y/SGx"^ + 20x + 17], 

— I < X < |. The point (x = |, y = — 1) is a classical tetracritical point at which the spiral 
phase V meets the FM phase, the striped collinear AFM phase, and the AFM Neel phase 
(and see Fig. 3 of Ref. [45| for further details). Thus, in our case, where J3 = J2 and Ji < 0, 
the classical spiral phase V exists in the narrow region < J2/\Ji\ < |- Naturally this 
includes the point J2/\Ji \ = y discussed above at which the FM and striped collinear AFM 
phases would meet in the absence of the spiral phase V as a stable gs phase intervening 
between them. 

In all of our results below for the FM honeycomb system we henceforth set 

Ji = — 1 with no loss of generality, since this simply sets the overall scale of the Hamiltonian, 
and we consider the case where J3 = J2 > 0, such that both the NNN and NNNN bonds 
act to frustrate the ferromagnetism. 

III. COUPLED CLUSTER METHOD 



The CCM (see, e.g., Refs. 60N62| and references cited therein) that we use here is one of 



the most powerful and most versatile modern techniques in quantum many-b ody theory. It 



has been used to study various quantum magnets (see, e.g., Refs. Ill 4l4lj6ll -l71| and refer- 
ences cited therein) very successfully. The method is particularly suitable for investigating 
frustrated systems, due to the fact that some of the main alternative methods are restricted 
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by certain problems that arise in such cases. For instance, quantum Monte Carlo (QMC) 
techniques suffer from the infamous and well-known "sign problem" for such systems. The 
exact ED method is also usually restricted by available computational power to relatively 
small finite-sized lattices. Nevertheless it can often be used, as here, to provide a handy tool 
to check and validate the results of other numerical or approximate methods. 

We briefly describe here some of the important features of the CCM as applied to spin- 



lattice problems (and see, e.g., Refs. Il|j61rl65| and references cited therein for further 



details). The starting point for any CCM calculation is to select a normalized state |$) 
as a reference or model state against which to incorporate in a systematic and potentially 
exact fashion the correlations present in the exact ground state. We often use a relevant 
classical ground state as the model state for spin systems for the sake of convenience, but 
other appropriate states may certainly also be used. In order to treat each site equivalently, 
a mathematical rotation of the local axes of the spins is conveniently performed in such a 
way that all spins in the reference state align along the same direction, say the negative 
z-axis. Clearly, such rotations leave unchanged the SU(2) commutation relations between 
components of the spin operators. 

The exact ket and bra gs energy eigenstates, |\&) and of the many-body system are 
then parametrized in the CCM form as: 



|vl/)=e^|$); S = J2SjCt, (3) 

(§| = (<l>|^e-^; S =l + J2^iC7, (4) 

where 

H\^)=E\<if); {^\H = E{^\, (5) 

are the Schrodinger gs ket and bra equations respectively. The multiconfigurational creation 
operators C/ = {CTV defined so that ($|C/ = = Cyl^) ; V/ 7^ 0, and where we have 
defined Cq = 1 = Cq . They are required to form a complete set of mutually commuting 
many-body creation operators in the Hilbert space, defined with respect to |$) as a cyclic 
vector. Clearly the states are normalized such that (^E'l^E') = (*^'|^) = ('^1'^') = 1- For 
spin-lattice systems they take the form of multi-spin raising operators and are written as 
products of single-spin raising operators, = s'j[s'^^ ■ ■ ■ s'^'^, where = Sj + ^sj- The 
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TABLE I: Number of fundamental LSUBm configurations (Nf) for the collinear striped 
AFM state of the spin-i Ji-J2^Jz honeycomb model. 

Method Nf 
LSUB6 72 
LSUB8 941 
LSUBIO 14679 
LSUB12 250891 



gs energy is calculated in terms of the correlation coefficients {Si} as E = (\E'|if |^) = 
($|e~'^ife"^|$); and the average on-site magnetization M in the rotated spin coordinates 
is calculated equivalently in terms of the coefficients {Si,Si} as M = — -^(^| J2j^=i^j\'^)y 
which now plays the role of the order parameter. Finally, the complete set of unknown 
ket- and bra-state correlation coefficients {5/, 5/} is calculated by requiring the expectation 
value H = (\['|iJ|\E') to be a minimum with respect to all parameters {Sj, 5/ ; V/ 7^ 0}. This 
readily leads to the coupled set of nonlinear equations for the ket-state creation correlation 
operators {Sj}, {^\CYe~^He^\^) = ; V/ 7^ 0, and to the coupled set of linear equations, 
{^\S{e~^ He^ — E)Cj'\^) = ; V/ 7^ 0, which can then be solved for the bra-state destruction 
correlation operators {Sj}. 

When all many-body configurations {/} are included in the expansions of the correlation 
expansions operators S and S, the CCM formalism is exact. However, it is necessary of 
course in practice to use approximation schemes to truncate the sets of configurations {/} 
contained in the expansions of Eqs. ([3]) and (j4]) for the CCM correlation operators. For 
systems defined on a regular periodic spatial lattice as here, it is convenient to use the well- 
established LSUBm approximation scheme in which all possible multi-spin-flip correlations 
over different locales on the (here, honeycomb) lattice defined by m or few contiguous lattice 
sites are retained. Clusters are defined to be contiguous in this sense if every site in the cluster 
is adjacent (as a nearest neighbor) to at least one other site in the cluster. This is the scheme 
we use for all our results presented below. The number A'^^ of independent fundamental 
clusters (i.e., those that are inequivalent under the symmetries of the Hamiltonian and of 
the model state) increases rapidly with the truncation index m, as shown in Table [T] for 
the present spin-| Ji-J2^J3 model on the honeycomb lattice, where we use the natural 
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lattice geometry itself to define the notion of adjacency inherent in the definition of the 
LSUBm scheme. We see, for example, that the number Nf of such fundamental clusters 
(and hence the number of simultaneous nonlinear equations we need to solve for the retained 
correlation coefficients {5/}) for the striped model state is 250891 at the highest LSUB12 
level of approximation that we utilize here. The corresponding numbers, A'^^, of fundamental 
configurations are appreciably higher at a given LSUBm level when the spiral phase V is 
used as the CCM model state, due to the considerably reduced symmetry. It is necessary 
to use massive parallelization and supercomputing resources in order to perform the CCM 
calculations at such high level of approximation.— Thus, for example, to obtain a single data 
point (i.e., for a given value of J3 = J2) for the striped model state at the LSUB12 level 
typically requires about 0.5 h computing time using 1000 processors simultaneously. 

We present CCM results below based on the striped coUinear AFM state as model state, 
at various LSUBm levels of approximation with m = {6,8,10,12}, and also in the cor- 
responding m -7- 00 extrapolated limits (LSUBoo) based on the well-tested extrapolation 
schemes described below and in more detail elsewhere.— "—'^'^ We have also performed ex- 
trapolations for the data set with m = {6,8, 10}. Both sets of results agree well with one 
another, which gives added credence to our results. Note that we do not use the LSUBm 
approximation scheme for values m < 6 of the truncation index, since these low-order ap- 
proximations will not capture the natural hexagonal structure of the lattice. We remark that, 
as always, the CCM exactly obeys the Goldstone linked-cluster theorem at every LSUBm 
level of approximation. Hence we work from the outset in the limit N ^ 00, where is the 
number of sites on the honeycomb lattice, and extensive quantities such as the gs energy are 
always guaranteed to be linearly proportional to in this limit. 

We clearly do not need to perform any finite-size scaling of our results, as all CCM 
approximations are automatically performed from the outset in the infinite- lattice limit, 
A^ —J- 00, as discussed above. It is, however, necessary to extrapolate to the exact m -> 00 
limit in the LSUBm truncation index m, in which limit the complete (infinite) Hilbert space 
is reached. For the gs energ y p er spin, E/N, a. well-tested and very accurate extrapolation 



ansatz (and see, e.g., Refs., |12lJl3lJ39lJ64lj73lJ74| ) is 



E (m) / N = tto + aim ^ + a2m ^, (6) 
while for the magnetic order parameter, M, different schemes have been employed for dif- 
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ferent situations. For models showing no or only relatively small amounts of frustration, a 



well-tested and accurate rule (and see, e.g., Refs. |64jj73| ) is 



M{m) = bo + feim"^ + b2m~^ . (7) 

For highly frustrated systems, particularly those showing a gs order-disorder transition, a 
more appropriate extrapolation rule with fixed exponents that has been found to give good 



results (and see, e.g., Refs. |l2lJ39| ) is 



M(m) = Co + + 62"^"^^^ • (8) 

We give illustrations here of the use of each of these schemes, wherever and whenever 
possible. 

IV. RESULTS AND DISCUSSION 

We now present and discuss our CCM results. In order to have an independent check 
on the accuracy and consistency of our CCM results, we have also performed additional 
computations of selected gs properties of the present models using the ED technique that is 



a well-establish 
e.g., Refs. [39 



led and successful tool for studying frustrated quantum spin systems (and see. 



48 



75 



78j). In Fig. [2]^a) we show the CCM results for the gs energy per spin, 
E/N, in various LSUBm approximations based on the striped state as CCM model state, 
as well as the exact gs energy for a finite lattice of size = 32. We also show separately, 
in Fig. Mjo), the extrapolated (m 00) results obtained from Eq. ([6]) using the data set 
m = {6, 8, 10, 12}. Comparison is made with the results for the corresponding AFM version 
of the model with Ji = +1. 

The CCM LSUBm data displayed in Fig. |2I^a) show that the gs energy results converge 
extremely rapidly as the truncation index m is increased, such that the difference between 
the LSUB12 results and the extrapolated (m — )■ 00) results obtained from Eq. is very 
small indeed. We note too that, just as in the corresponding AFM case of the model 
with Ji = +1, the various CCM LSUBm solutions based on the striped model state now 
also terminate at some lower termination point J2 as J2 is decreased. Such terminations 
of CCM solutions are very common and have been very well documented.— In all such 
cases a termination point always arises due to the solution of the CCM equations becoming 
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(a) m = 6, 8, 10, 12 (b) m cx) 

FIG. 2: (Color online) Ground-state energy per spin, E/N, versus J2 for the striped phase 
of the spin-i Ji-J2^Jz honeycomb model with Ji = — 1 compared with those for Ji = +1 

(with J3 = J2 > for both cases). The CCM results using the striped model state for 
various LSUBm approximations with m = {6, 8, 10, 12}, together with the ED (A^ = 32) 
results are shown in (a). We show in (b) the extrapolated CCM LSUBoo results using the 

m = {6, 8, 10, 12} data points fitted to Eq. ([6]). In all cases curves without symbols 
attached refer to the case Ji = —1, whereas the corresponding curves with symbols refer to 
the case Ji = +1. The FM result from Eq. ([2]) with s = | is also displayed. 



complex at this point, beyond which there exist two branches of entirely unphysical complex 
conjugate solutions.— In the region where the solution reflecting the true physical solution 
is real there actually also exists another (unstable) real solution. However, only the shown 
branch of these two solutions reflects the true (stable) physical ground state, whereas the 
other branch does not. The physical branch is usually easily identified in practice as the one 
which becomes exact in some known (e.g., perturbative) limit. This physical branch then 
meets the corresponding unphysical branch at some termination point with infinite slope, 
beyond which no real solutions exist. The LSUBm termination points are themselves also 
reflections of the quantum phase transitions in the real system and may be used to estimate 
the position of the phase boundary,— although we do not do so for this critical point in the 
FM model, since we have more accurate criteria that we now discuss. 

We note first from Fig. [2](a) that the LSUBm termination points using the striped state 
as the CCM model state for the present FM version of the model with Ji = —1, lie very 
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close indeed to the points where the curves cross (or nearly cross) the corresponding curve 
for the FM state given by Eq. ([2]). This gives us our first evidence that either there is 
no intermediate phase between the quantum striped phase and the FM phase for the case 
Ji = —1, or, if one exists, it can occur only over a very narrow regime indeed. This 
situation may be contrasted with that of the AFM version of the model (Ji = +1),— 
where the LSUBm results for the gs energy using the striped model state terminate before 
they meet the corresponding results using the Neel state as model state (which themselves 
also terminated at some upper termination points that were lower in value than the lower 
termination points for the striped state). In the latter case there is an intermediate plaquette 
valence-bond crystal (PVBC) phase. 

At the classical level the difference in the values of the gs energy per spin of the coUinear 
striped states between the two s = | cases (i.e., for positive and negative values of Ji with 
I Ji| = 1) is 0.25, independent of J2 and J3. The quantum versions follow this pattern for 
larger values of J3 = J2, as seen from Fig. |2l but the constancy in the difference breaks 
down at around J2 ~ 0.6, where the AFM case (Ji = +1) exhibits a critical point marking a 
transition to the PVBC phase, which then in turn undergoes a further phase transition to the 
Neel phase at another lower critical value. The corresponding best available COM estimates 
for those two critical values for the AFM case of the model with Ji = +1 are ~ 0.60 
and ~ 0.47 respectively.— In the present FM case of the model with Ji = — 1 we see 
no evidence (apart from the seeming termination of the solutions to the equations for the 
LSUB12 approximation based on the striped state as CCM model state very slightly before 
the gs energy crossing point with the FM state) of any similar intermediate state between 
the FM state and the collinear striped AFM state. If any such intermediate state exists at 
all, however, it must be confined to a very narrow region indeed around J2 ~ 0.11, probably 
confined to 0.10 < J2 ^ 0.12. We return to a more detailed discussion of this region later. 
For the moment we note only that it is much reduced from the region 0.1 < J2 < 0.2 in 
which the corresponding classical version of the model has the spiral phase V as its stable 
gs phase. 

For the present FM case with Ji = —1 the CCM LSUBm gs energy curves using the 
striped model state cross the corresponding gs energy curve for the FM state from Eq. ([2]) 
for m = {6, 8, 10} at corresponding critical values J2(LSUB6) ~ 0.1106 (where E^^jj^q/N ^ 
-0.2506), J2"(LSUB8) ^ 0.1101 (where ^Esubs/^ ~ -0.2511), and J|(LSUB10) ^ 0.1098 
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(where -E'lsubio/^ ~ —0.2515). The corresponding LSUB12 result for the gs energy using 
the striped state as the CCM model state appears to terminate just before meeting the 
gs energy curve for the FM phase. However we note that for such very high-order CCM 
calculations it becomes very computationally expensive to determine the termination point 
with high accuracy. If we use the extrapolated LSUBoo results for the gs energy for the 
striped phase by making use of Eq. and employing the whole data set m = {6, 8, 10, 12}, 
we thus need to perform a further very small extrapolation of the CCM results to lower values 
of J2 to find the presumed crossing point of the energies of the striped and FM phases, in 
the scenario in which these two phases meet at a first-order transition with no intermediate 
phase (that would itself be confined to the very narrow intervening region 0.10 ^ J2 ^ 0.12, 
as discussed above). As expected, simple power-law expansions give very accurate fits, and 
give crossing points very close to those above. Putting all the energy data together, our 
best estimate for the critical point of the first-order phase transition from the coUinear 
striped phase to the FM phase (in the scenario where this transition occurs directly, with no 
intermediate phase confined to the narrow region 0.10 < t/2 ^ 0.12) for the spin-| J1-J2-J3 
Heisenberg ferromagnet (with Ji = —1) on the honeycomb lattice, and with J3 = J2 > 0, is 
J| = 0.1095 ± 0.0005, at which point the gs energy per spin is E^N = -0.2518 ± 0.0006. 

We see from Fig. [2]^a) that the agreement between the ED (A^ = 32) and the CCM 
energies is very satisfactory. Moreover, due to the finite-size scaling of the gs energy, E/N = 
Co — a/N'^/'^ with a > (and see, e.g. Refs. [3] and [t?]), the difference between the CCM 
and the ED gs energies would become even smaller if finite lattices of larger size could be 
considered. The ED turnover point in the energy curve that marks the termination of the 
FM phase occurs at a value of about 0.1003 for the = 32 lattice used, and for the same 
reasons as above this value will increase as A^ is increased. Thus, in summary, while the 
CCM estimates for the gs energy per spin for the spin-| Ji-J2^Jz Heisenberg model on the 
extended infinite honeycomb lattice are much more accurate than the ED results, the latter 
do serve as an independent check on the former. 

The hypothetical phase transition (i.e., when the existence of the intervening spiral phase 
V is momentarily ignored) from FM order to collinear striped AFM order for the classical 
version of the FM model with Ji = —1 occurs at a value = 1 ^ 0.143, compared with 
the corresponding value J| ~ 0.110 found here. Thus quantum fluctuations act to stabilize 
the collinear AFM order, at the expense of the FM order, to higher values of frustration 
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than in the classical case. It is interesting to note that a similar situation was found in the 
FM version (Ji = — 1, J2 > 0) of the spin-| J1-J2 niodel on the square lattice,— where a 
quantum critical point exists at a value J| ~ 0.39 for a similar transition from FM order to 
collinear striped order, compared with a corresponding classical value of = 0.5. It is well 
known, from many cases studied, that quantum fluctuations almost always favor collinear 
states over noncoUinear ones (e.g., spiral or canted states). What is interesting in both the 
present case and the spin-i J1-J2 model on the square lattice cited above, is that quantum 
fluctuations seem also to favor one collinear state (namely the collinear striped AFM state 
in these two cases) where the quantum fluctuations are present, over another collinear state 
(namely the FM state in these two cases) where quantum fluctuations are absent. It is 
intriguing to wonder whether these are examples of a more general rule. 

We present results in Fig. |3] for the CCM collinear stripe order parameter M, as defined 
in Sec. IIIII Figure |3]^a) shows LSUBm results with m = {6, 8, 10, 12}, while Fig. [3]^b) shows 
the corresponding extrapolated CCM LSUBoo (m — )■ 00) results using both Eqs. ([7]) and 
(IH]). We note firstly that the CCM LSUBm order parameter results depend on the approxi- 
mation level m much more strongly than those for the gs energy. It is clear that the order 
parameter behaves similarly for large values of J2 for both the FM model (Ji < 0) and the 
AFM model (Ji > 0). However, once again there are considerable differences in the behavior 
of M between the two models for values of the frustration parameter J2/\Ji\ ^ 0.6. The 
extrapolated CCM results for M for the AFM model in Fig. [3]^b) clearly show the break- 
down of the quasiclassical collinear magnetic LRO near the critical value of <Z2^~ 0.6, i.e., 



significantly above the classical transition point = 0.5 (and see, e.g., Refs. 44j48|j49|j55|). 



Indeed, the CCM estimate for the critical value of the frustration parameter in the AFM 
case for the disappearance of collinear striped order is ~ 0.60 from the point at which M 
becomes zero, using the extrapolation scheme of Eq. (|8]).— By contrast, the order parameter 
for the FM model stays almost constant over the whole parameter region shown in Fig. [31 
We do not observe any indication of the breakdown of the collinear striped magnetic LRO 
until J2 ~ 0.11 for the FM model, which is below the hypothetical classical transition point 
J| fti 0.143, as we observed previously in the results for the gs energy. 

Lastly, we present results for various spin-spin correlation functions for the FM as well 
as for the AFM model in Fig. H Figure SJ^a) shows the CCM LSUBIO resuhs and Fig.H^b) 
shows the corresponding ED results. Once again we note that for large values of the frus- 
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FIG. 3: (Color online) Ground-state magnetic order parameter, M, for the striped AFM 



Ji = +1 (with J3 = J2 > for both cases). The CCM results for various LSUBm 
approximations with m = {6, 8, 10, 12} are shown in (a). We also show in (b) the 
extrapolated LSUBoo results using the m = {6, 8, 10, 12} data points. The curves labelled 
LSUBoo(l) and LSUBoo(2) are obtained using Eq. ([7]) and Eq. ([8]) respectively. In all 
cases curves without symbols attached refer to the case Ji = —1, whereas the 
corresponding curves with symbols refer to the case Ji = +1. 

tration parameter J2 the corresponding spin-spin correlations functions for both the FM 
(Ji = —1) and AFM (Ji = -|-1) models agree remarkably well with one another for both 
the CCM and ED calculations. Furthermore, for the FM model the agreement of the CCM 
correlation functions with the ED data is excellent. For the AFM model the agreement 
between the CCM and ED results is again excellent for values of J2 above the transition 
point at which the AFM coUinear striped order disappears, namely ~ 0.60, but around 
and below this value there are noticeable differences. In particular, the very steep change in 
the correlation functions at J2 ~ 0.62 present in the ED (A^ = 32) data for the AFM model 
is not observed in the CCM LSUBIO data. Instead the CCM data show a smoother change 
in that region. However, we have argued^ that for J < J^^ 0.60 no striped magnetic 
LRO order exists. Indeed we argued that no magnetic LRO order exists at all for the AFM 
model in the regime < J < J2^, where instead we have a PVBC state. Hence, it is 
not surprising that the CCM solution in a finite order of LSUBm approximation based on 



state of the spin-^ Ji-J2^J?, honeycomb model with Ji 



1, compared with those for 
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FIG. 4: (Color online) Selected spin-spin correlation functions (s(0) ■ s(r)) for the spin-| 

J1-J2-J3 model on the honeycomb lattice with J3 = J2, for both the FM case (with 
Ji = —1) and the AFM case (with Ji = +1), using (a) the CCM with the striped coUinear 
state as model state, at the LSUBIO level of approximation; and (b) the ED method on a 
lattice of size = 32. Values r = 1, 1.732 and 2 correspond respectively to NN, NNN, and 
NNNN pairs of spins. In all cases curves without symbols attached refer to the case 
Ji = —1, whereas the corresponding curves with symbols refer to the case Ji = +1. 



the collinear stripe reference state does not provide such accurate results for the correlation 
functions inside this magnetically disordered phase. 

To conclude, we return to examine more closely the very narrow region 0.10 ^ J2 ^ 0.12 
for which our CCM results based on the striped AFM state as model state could not exclude 
the possibility of an intervening phase between the striped AFM and the FM phases. In Fig. [5] 
we show a more detailed view of the ED results for the same spin-spin correlation functions 
shown in Fig. 111(b) in this narrow region just above the FM transition point. The ED data 
does definitely indicate the existence of a phase in precisely the region 0.10 < J2 ^ 0.12. 
It is difficult from this data to say with any certainty whether or not the state is the 
quantum-mechanical remnant of the classical spiral phase V that exists in the classical 
regime 0.1 < J2 < 0.2. Furthermore, without ED calculations on larger lattices, for which 
the computational cost would be prohibitive, it is also not possible to say whether these 
results over such a narrow region are an artefact of the finite lattice size. Our results are 
summarized in Sec. |Vl 
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FIG. 5: (Color online) An expanded view near the FM boundary of the same spin-spin 
correlation functions (s(0) ■ s(r)) shown in Fig. HJ^b) for the spin-| Ji-J2^^3 model on the 
honeycomb lattice for the FM case with Ji = — 1 and J3 = J2 using the ED method on a 

lattice of size = 32. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have presented results on the gs properties of the spin-| Ji-J2^Jz 
Heisenberg model with FM NN (Ji = —1) exchange bonds in the presence of frustrating 
AFM NNN (J2 > 0) and NNNN (J3 > 0) exchange bonds of equal strength (J3 = J2) on 
the honeycomb lattice, using both the CCM and Lanczos ED. By comparison with previous 
studies for the AFM (Ji = +1) version of the model,— we find similar behavior for both 
models for values J2 > 0.6|Ji|, but for values of J2 ^ 0.6|Ji| the models differ markedly. 
The results of the present paper for the FM version of the model and that of the previous 
paper— for the AFM version may conveniently be combined and summarised in the phase 
diagram shown in Fig. |6l 

We note that, by contrast with the corresponding model with AFM NN exchange (Ji = 
+1) we do not find indications for a non-classical magnetically-disordered phase for the 
model with FM NN exchange ( Ji = — 1). If such a phase exists at all it must be confined to 
a very small range of the frustration parameter around 0.10 ^ J2 ^ 0.12. However, any such 
phase is much more likely to be a quasiclassical remnant of the spiral phase V that exists in 
the corresponding classical model (with Ji = —1) in the parameter regime 0.1 < J2 < 0.2. 
As expected, quantum fluctuations then usually favor a collinear phase over a noncollinear 
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FIG. 6: (Color online) The phase diagram of the spin-i Ji-J2^^3 honeycomb model in the 
J1-J2 plane, for the case J3 = J2. The continuous transition between the AFM Neel and 
PVBC phases at J2I J\ ~ 0.47 is shown by a broken line, while the first-order transition 
between the PVBC and AFM striped phases at J2/J1 ~ 0.60 is shown by a solid line. Our 
results indicate that the transition between the striped AFM and FM phases is either a 
first-order one at J2/J1 ~ —0.11 or occurs via an intermediate phase, probably with 
noncollinear spiral order, which exists in the region —0.12 < J2/J1 ^ 0.10. The region 
between the FM and AFM Neel phases with J3 = J2 < has not been investigated by us. 

one, and the extent of any spiral phase is smaller in the quantum spin-^ case than in the 
classical (s — 0) case. 

In one scenario the results presented here for the case Ji = — 1 indicate a direct first- 
order transition between the two magnetically ordered phases, namely the FM ground state 
at small values of the frustration parameter J2 and the striped collinear AFM ground state 
at larger values of J2. Our best estimate of the phase transition point is then = 0.1095 ± 
0.0005. Although in this scenario a quasiclassical gs phase (viz., the collinear striped AFM 
state) exists in the whole parameter space down to the FM gs phase, the frustration might 
still have a strong effect on the low-temperature thermodynamics near the transition point 
at J| = 0.1095 ± 0.0005.-^'^ For values J2 > J2 the FM multiplet becomes a low-lying 
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excitation, and this might lead to an additional low-temperature peak in the specific heat 
C(T).— 1^1^ We note that indications for such an additional low-temperature peak in C{T) 
were also found on the FM side near such a transitional^ (i.e., at J2 ^ </|) in other frustrated 
spin models. 

In an alternative scenario our results also indicate the possibility of an intervening phase 
between the coUinear FM and striped AFM phases. Any such phase, however, is limited to 
lie within the very narrow range 0.10 ^ J2 ^ 0.12, as shown in Fig. [6l In principle we could 
more accurately establish the existence of such a phase as a quasiclassical remnant of the 
classical spiral phase V, and thence also more accurately establish its phase boundaries, by 
performing another comparable set of CCM calculations to those performed here with the 
striped AFM state as model state, but using instead the spiral state V as model state. Such 
calculations would be much more onerous and computationally expensive, however, since 
on the one hand the number Nf of fundamental CCM configurations at a given LSUBm 
level is greater for the spiral model state than for the striped model state and, furthermore, 
the CCM results would need to be optimized at a given LSUBm level with respect to the 
spiral pitch angle parameters by minimizing the corresponding result for the energy per spin 
separately for each set of values for the bond strength parameters. 

We note finally that we have not yet investigated the present model in the case where 
J3 = J2 < 0. For the FM version of the model when Ji < also, the FM phase is then 
obviously the stable ground state. Conversely, when Ji > and frustration occurs, there 
is a direct first-order transition in the classical version of the model between the FM and 
Neel AFM states at a value J2/J1 = —1- Following the discussion in Sec. IIVI we might 
expect that quantum fluctuations could again act either (a) to retain the direct transition 
but to stabilize the coUinear AFM order in preference to the FM order, thus pushing the 
phase boundary to a somewhat lower value, J2/ Ji < —1, for the spin-| case; or (b) to 
permit an intervening state with no classical counterpart. Indeed, very preliminary CCM 
calculations indicate that scenario (a) is realized and that this corresponding critical point 
may be pushed to a value J2/J1 ~ —1.15 ± 0.05. We hope to report in more detail on this 
region and to give a more accurate value of this phase boundary in a future paper. 

As discussed briefly in Sec. HI it has been proposed^'^*^'^ that the competition between 
FM Heisenberg interactions between NN pairs of spins and AFM interactions between other 
spins in frustrated spin-| systems on the square lattice could lead to gapless spin-liquid 
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states with multipolar order (e.g., spin-nematic states) adjacent to the FM state. Similar 
states have also been proposed to arise in frustrated multiple cyclic spin-exchange models 
on the triangular lattice with FM NN pairwise interactions,— either in the presence of a 
magnetic field (where octupolar order occurs) or in its absence (where quadratic or nematic 
ordering occurs in a state bordering the FM state). In the case of the frustrated honeycomb- 
lattice ferromagnet considered here we have found no evidence for such states. However, 
the multipolar-ordering phenomenon in the zero-field case considered here is evidently rather 
fragile, and in the square-lattice case for the spin-i FM version of the J1-J2 model (i.e., with 
Ji < 0) even their existence has been questioned in recent rather accurate work— that also 
employed both high-order CCM and ED techniques. No evidence was found for such states 
either in a very recent Schwinger boson study on the square lattice,— using the same FM 
version of the spin-^ Ji-J2^^3 Heisenberg model that we studied here on ther honeycomb 
lattice. Nevertheless, the history of the study of quantum magnets has shown us that the 
detection of phases with novel quantum ordering, such as nematic states of various kinds, 
is extremely subtle. In particular, the present honeycomb-lattice model surely warrants 
further investigation before the absence of nematic states in the FM case discussed here is 
considered definite. 

Finally we mention that frustrated ferromagnets are also interesting with respec t to multi- 

Sal-M). The 
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magnon bound states appearing in high magnetic fields (and see, e.g., Refs. 
present model also warrants further investigation when the coupling to an external magnetic 
field is included. 
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